Topologically protected quantum transport in locally exfoliated bismuth at room 

temperature 
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We report electrical conductance measurements of Bi nanocontacts created by repeated tip-surface 
indentation using a scanning tunneling microscope at temperatures of 4 K and 300 K. As a function 
of the elongation of the nanocontact we measure robust, tens of nanometers long plateaus of conduc- 
tance Go = 2e^ /h at room temperature. This observation can be accounted for by the mechanical 
exfoliation of a Bi(lll) bilayer, a predicted QSH insulator, in the retracing process following a tip- 
surface contact. The formation of the bilayer is further supported by the additional observation of 
conductance steps below Go before break-up at both temperatures. Our finding provides the first 
experimental evidence of the possibility of mechanical exfoliation of Bi bilayers, of the existence of 
the QSH phase in a two-dimensional crystal, and, most importantly, of the observation of the QSH 
phase at room temperature. 



Topological insulators present a gap in the bulk, but 
host surface states protected against backscattering by 
time reversal symmetry [l|. This implies that they are im- 
mune to non- magnetic disorder-induced localization, i.e., 
they are able to carry electrical current on the surface re- 
gardless of imperfections. 2D TIs[2], actually predicted 
before their three-dimensional (3D) counterparts [35, are 
expected to exhibit the so-called quantum spin Hall 
(QSH) phase, a spin filtered version of the integer quan- 
tum Hall effect [J] . While the most exotic experimental 
manifestation of this phase is through a nearly universal 
spin Hall conductivity of ~ e/27r, a truly universal charge 
transport is expected to manifest, e.g., as a two-terminal 
conductance Gq = 2e^ /h. 

To date, two types of 2D systems have been predicted 
to be QSH insulators: two-dimensional crystals such as 
graphene[2j or Bi(lll) bilayers [5j and semiconductor het- 
erojunctions such as CdTe/HgTe[6] or, more recently, 
InAs/GaSb quantum wells [7]. Transport measurements 
in CdTe/HgTe[6] and InAs/GaSb[8] quantum wells have 
revealed the presence of protected edge states and pro- 
vided the first experimental evidence of the QSH phase to 
date. The QSH state in 2D crystals, on the other hand, 
has not been experimentally confirmed to date. The fact 
that spin-orbit coupHng (SOC) in graphene is so weak 
precludes the observation of the QSH phase in this ma- 
terial. Bismuth, on the contrary, presents a naturally 
large SOC, its mechanical and electronic properties are 
well characterized for bulk and surface [9], in nanowire 
form[lO|, and, recently, the existence of edge states in 
Bi(lll) bilayers has been reported [11]. Other proposals 
stay, at this moment, at a more speculative level[12]. 

Cleavage techniques are becoming common in the 
quest for 2D crystals [l3|. With a few exceptions [14], they 
remain largely unexplored in the field of topological in- 
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FIG. 1. (a) Selected curves obtained during the process of 
rupture of a Bi nanocontact formed after indentation of a 
tip on a surface at 4 K. (b) Detail of a selected curve where 
four conductance sub-plateaus can be appreciated below Go- 
The insets represent possible atomic configurations compati- 
ble with the sub-plateau values. 



sulators (TI's). We report here, using scanning tunnel- 
ing microscope (STM) based mechanical and electrical 
characterization techniques, the first evidence of the QSH 
phase in a two-dimensional crystal such as an exfoliated 
Bi(lll) bilayer. The procedure follows closely that of 
creation of atomic size metallic contacts[15], namely, re- 
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peated indentations of a tip against a surface while mea- 
suring the current at low voltages (typically 100 meV). 
The resulting conductance traces contain information on 
the atomic and electronic structure of the contact. The 
experiments where done at room temperature (and am- 
bient conditions) and at low temperatures (4K) in cryo- 
genic vacuum using two Bi samples of different origin [ 
Related experiments have been reported in the past[ 
ll9[, lacking a convincing interpretation of the results, but 
have, in part, motivated our work. 

As a rule of thumb (in particular for monovalent el- 
ements) each atom forming the minimum cross section 
of a metallic nanocontact contributes to the conductance 
with a value in the vicinity of Go[l5[. Thus, on retracing 
the tip after indentation, one-atom cross-section contacts 
formed right before complete break-up are signaled by 
a conductance plateau (as a function of the STM piezo 
elongation) near the quantum of conductance Gq. Unlike 
their more common metallic counterparts, the conduc- 
tance traces of Bi nanocontacts at 4 K present several 
sub-plateaus below Go before break-up (see Fig. la). 
In Fig. lb, we have singled out one trace that mani- 
festly shows four sub-plateaus before the final rupture. 
Since each plateau corresponds to an elastic deformation 
of an atomic configuration in between plastic deforma- 
tion events, one can only conclude that few-atom cross 
section Bi nanocontacts (as the ones shown in the in- 
sets) are much less conducting than their metallic coun- 
terparts. Being a common practice to record conductance 
histograms to statistically characterize nanocontacts |15J. 
we have done so using thousands of rupture traces. We 
find a large statistical weight below Go, but no charac- 
teristic values can be appreciated [16|. 

In air at 300 K the mechanical stability of the exper- 
imental set-up is inferior to that at low temperatures, 
making it more difficult to carry out measurements in 
a systematic way. Previously, the samples were fully 
characterized[16J, paying attention to surface contami- 
nation and, in particular, to oxidation. This turned out 
to be very small, as long known for Bi[20]. Starting from 
conductance values at a maximum indentation of the or- 
der of lOGo, most of the traces either do not break for 
the allowed range of the piezo, typically up to 300 nm, 
or fall to zero showing no apparent reproducible behav- 
ior. Starting from smaller indentations (up to 3-6 Go), 
however, consecutive traces may repeatedly show small 
features or short plateaus, including plateaus at Go (see 
Fig. 2a) which were completely absent in the 4 K mea- 
surements. Notice also the different length scales in Figs. 
la and lb (1 nm) and Figs. 2a and 2b (100 nm). Most 
notably, long conductance plateaus appear on occasions 
on some traces (see Fig. 2b), the length of these reach- 
ing up to hundreds of nanometers. While some of them 
correspond to the initial elongation elastic process (the 
ones around 3Go), most of them appear after a few plas- 
tic events. One way to quantify the information in these 




50 100 150 

Piezo displacement (nm) 



o 
0) 

o 

c 

5 2 
o 

1 1 

O 

o 



a 





K. 


t^-- 


5 


t5 







1000 



100 200 300 

Piezo displacement (nm) 



(/) 



.9 500 




Conductance (G ) 



FIG. 2. (a) Examples of conductance recorded during the 
process of rupture of a Bi nanocontact formed by repeated 
indentation of a Bi tip on a Bi surface at 300 K. (b) Conduc- 
tance traces exhibiting plateaus with elongations of tens and 
hundreds of nanometers attributed to the occasional exfolia- 
tion of large flakes of Bi(lll) bilayers between tip and sub- 
strate. Most of the plateaus are pinned at Go- (c) Histogram 
of conductance obtained from tens of traces exhibiting long 
plateaus. The histogram exhibits a clear peak at Go which 
reveals that this value is particularly frequent and robust. 



traces is to plot a conductance histogram only with them 
(we have selected those exhibiting plateaus longer that 
50 nm regardless of the plateaus conductance values). 
The result is incontrovertible, showing a large peak at 
Go since most of the plateaus are largely pinned at the 
quantum of conductance (see Fig. 2c). 

An exact conductance value of Go is generally expected 
to occur when, first, the Fermi wavelength is compara- 
ble to the constriction width and, second, the constric- 
tion potential is adiabatic in the current direction. This 
happens, for instance, in a gated two-dimensional elec- 
tron gas[21| or in monovalent metal atomic contacts |15|. 



Crystalline bulk Bi presents two types of carriers (elec- 
trons and holes) with long Fermi wavelengths (tens of 
nm) and unusually small masses [22|. A first scenario, al- 
ready put forward to explain the early experiments in 
Ref. 18, invokes the formation of large cross section 
three dimensional constrictions (thousands of atoms) so 
that the lateral movement of the long-wavelength elec- 
trons gets quantized. It is, however, difficult to predict 
how the bulk electronic structure can carry over to con- 
strictions presumably lacking long crystalline order. In 
addition, adiabaticity imposes severe restrictions to the 
shape of the constrictions. Furthermore, it is quite hard 
to imagine how these constrictions may be stable during 
very long piezo displacements. Even more difficult to ac- 
cept is the fact that they must exist right before rupture 
without then going through a very large number of plastic 
deformation events. A second scenario where atom-size 
constrictions are responsible for the appearance of Go is 
ruled out by the observation of sub-quantum plateaus in 
most traces (see Figs. 1 and 2) and also by our calcula- 
tions (see below). 

We now give our rationale to our findings. Bulk Bi can 
be viewed as a layered material in the (111) direction ^16|; 
much as graphite, but with a stronger electronic cou- 
pling between layers. In fact, a Bi(lll) bilayer is similar 
to graphene in many respects, in particular the atomic 
structure is that of a puckered honeycomb lattice, with 
the atoms of one sublattice shifted perpendicularly with 
respect to those of the other sublattice. A Bi(lll) bilayer 
has been predicted to be a QSH insulator [5| with a large 
bulk gap and three pairs of helical edge states [23]. This 
odd number of edge states guarantees that the transport 
properties of a Bi(lll) flake should be robust against 
size, shape, and weak disorder. In particular, it should 
generically exhibit a universal two-terminal conductance 
of Go. 

Given the layered structure of Bi, the general picture 
that naturally emerges from our experiments is that of a 
contact-induced exfoliation of a Bi(lll) bilayer (see Figs. 
3a and 3b). We first note that the inter-bilayer coupling 
is about ten times weaker than the intra-bilayer cova- 
lent bonding [9|. The Bi(lll) surface presents a Debye 
temperature in the 70-80 K range |24j so that at ambient 
conditions it should not come as a surprise that, once 
a gentle contact is made on the appropriate surface ori- 
entation, retracing the tip can peel off a single bilayer. 
One extreme of this bilayer (while the rest remains in 
contact with the surface) sticks to the tip which now 
acts as an electrode (see Figs. 3a, 3b and Supplemen- 
tal Material movie [16,]). The specifics of the proposed 
mechanical exfoliation are unknown to us, being unclear 
whether the tip breaks the surface or simply elastically 
deforms it after indentation, and whether slip directions 
on the (111) planes favor slipping over detaching mech- 
anisms. The number of atoms participating in the con- 
tact at maximum conductance is also unknown, but it 
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FIG. 3. (a) Pictorial representation of the proposed process 
of exfoliation of a Bi(lll) bilayer (in red) after contact with 
the STM tip. One (out of three) helical edge channel is also 
represented, (b) Same as in (a) but for the last stages of 
piezo elongation before break-up where a few-atom constric- 
tion forms, (c) Band structure of Bi(lll) bilayer nanorib- 
bons along ZZ (left) and AC (right) directions as obtained 
from a self- consistent tight-binding model. The insets show 
a top-view section of these nanoribbons. The edge bands 
connecting valence and conduction bands are shown in blue. 
The Fermi energy (set to zero) crosses three times these edge 
bands which gives rise to three conduction channels. 



remains large enough after exfoliation so that the theo- 
retically expected intrinsic conductance Go of the bilayer 
is not masked by the contact resistance. 

With this picture in mind, the conductance sub- 
plateaus below Go can only reflect the breaking of the 
tip-bilayer contact or the progressive breaking of the bi- 



layer, a process in which a bidimensional nanocontact 
may form (see Fig. 3b and insets in Fig. 4a). Since a dis- 
ordered Bi(lll) bilayer is expected to exhibit a maximum 
conductance of Go, we can only expect the nanocontacts 
thus formed to conduct less due to inter-edge backscat- 
tering at the constriction (this is more clearly observed 
at 4 K). Noticeably, at 4 K, long plateaus do not appear 
near Go . We attribute this absence to the stronger inter- 
bilayer binding at low temperatures which prevents the 
exfoliation of a sufficiently large flake. Interestingly, the 
up-turns of the sub-plateaus rarely surpass Go- 

We now support our hypothesis with atomistic cal- 
culations. We model Bi with a 4-orbital tight-binding 
Hamiltonian with the parametrization given by Liu and 
Allen [22]. In order to account properly for the electronic 
structure at the edges, we carry out self-consistent calcu- 
lations including local electron-electron interactions [16|. 
As can be seen in Fig. 3c, Bi(lll) bilayer nanoribbons 
grown along zigzag (ZZ) and armchair (AC) crystallo- 
graphic orientations present a gap with valence and con- 
duction bands connected in a continuous manner by an 
edge state band. This edge band cuts three times the 
Fermi energy and presents the same topology as that 
computed with more sophisticated DFT calculations [23|, 
which gives us confidence in our model. We have also ver- 
ified that the nature of the wave functions in this band 
follows that expected in a QSH insulator |16|. 

Equipped with a reliable atomistic description of Bi, we 
now perform coherent quantum transport calculations. 
(While inelastic effects cannot be entirely ruled out, bulk 
Bi carriers are known to have very long mean free paths 
even at room temperature [25].) We introduce constric- 
tions in the nanoribbons emulating the breaking of the 
bilayer flake (see insets in Fig. 4a). For simplicity's sake 
we have chosen to perform this type of calculation in- 
stead of the one for a contact between a tip (of unknown 
shape) and the flake (as shown in Figs. 3a and 3b). A 
common way to calculate the conductance (as a function 
of energy) G{E) is through the Green's function formal- 
ism and the partitioning technique[l6, 26, 27], now tak- 
ing fully into account the atomic SOC^, |28i]. Results of 
G{E) for four different constriction geometries in AC and 
ZZ nanoribbons are shown in Fig. 4a. Figure 4b shows G 
at the Fermi energy {Ep = 0) for these four geometries 
as a function of the number of bonds at the constriction 
(i.e., as a function of the width). Only in disorder- free 
ribbons with perfectly defined edges the three channels 
can fully transmit. For constriction widths just two or 
three bonds smaller than that of the ribbons the conduc- 
tance drops in all cases to values in the close vicinity of 
Go- Only constrictions featuring less than approximately 
five bonds conduct below Go- (The robust appearance of 
Go also holds for other types of disorder [16.].) These re- 
sults can be simply understood in terms of strong intra- 
edge backscattering of two of the edge channels and ab- 
sent or partial inter-edge backscattering of the remaining 
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FIG. 4. (a) Conductance vs. electron energy for four con- 
striction models in ZZ and AC Bi(lll) bilayer nanoribbons. 
(b) Conductance at the Fermi energy for many cases as those 
in (a) as a function of the number of bonds remaining in the 
constriction. The sub indexes denote central- (c), lateral- (1), 
and wedge- type (w) constrictions. 



helical pair at the narrowest section of the constriction 
(see Figs. 3a and 3b). Our estimate of a minimum of 
five bonds for a constriction to support G{Ef) = Go is 
consistent with the maximum number of plastic events 
typically seen in the experiments below Go before the 
constriction definitely breaks (see, e.g.. Fig. lb). 

In summary, we have reported and offered a consistent 
interpretation to the appearance of extremely long quan- 
tum of conductance plateaus in the breaking process of 
Bi nanocontacts. We attribute it to the occasional lo- 
cal exfoliation of Bi bilayers, predicted to be 2D topo- 
logical insulators. Other interpretations cannot be ruled 
out, but should account for three extraordinary facts: (i) 
the quantum of conductance cannot be associated with 
a single- atom contact, but with a nanoscopic constric- 
tion, (ii) it appears at room temperature, and (iii) it is 
robust to small concentration of contaminants since the 
experiments are performed in air. 



The first two authors, CS and DGM, contributed 
equally to this work. The former on the experimen- 
tal side and the latter on the theoretical part. We 
thank the Spanish MICINN for financial support through 
Grant Nos. FIS2010-21883, FIS2011-23488, and CON- 
SOLIDER CSD2007-0010 and Generalitat Valenciana 
through Grant Nos. ACCOMP/2012/127 and PROM- 
ETEO/2012/011. DG acknowledges the Centro de Com- 
putacion Cientifica at UAM for computational support. 
Cristina Almansa is gratefully acknowledged for her help 
with the transmission electron microscope. 



[1] 

[2] 

[3] 

[4] 

[5] 
[6] 

[7] 

[8] 

[9] 
[10] 



M. Z. Hasan and C. L. Kane, [Rev7 Mod. Physl 82, 3045 
(2010). 

C. L. Kane and E. J. Mele, Mi^. Rev. Lett. 1 95, 226801 
(2005). 

H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. 
Zhang, Nature Physics 5, 438 (2009). 
The Quantum Hall Effect, edited by R. E. Prange, S. M. 
Girvin, and K. v. Klitzing (Springer, 1986). 
S. Murakami, Phys. Rev. Lett. 97, 236805 (2006). 
M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, 
ISciencel 318. 766 (2007). 

C. Liu T. L. Hughes, X.-L. Qi, K. Wang, and S.-C. 
Zhang, [Phys. Rev. Lett.] 100, 236601 (2 008). 
L Knez, R.-R. Du, and G. Sullivan, P hys. Rev. Lett.| 107, 
136603 (2011). 

P. Hofmann, Progress in Surface Science 81, 191 (2006). 
Z. Zhang, D. Gekhtman, M. S. Dresselhaus, and J. Y. 
Ying, IChem istry of Materials 11, 1659 (1999). 
[11] F. Yang, L. Miao, Z. F. Wang, M.-Y. Yao, F. Zhu, Y. R. 
Song, M.-X. Wang, J.-P Xu, A. V. Fedorov, Z. Sun, G. B. 
Zhang, C. Liu, F. Liu, D. Qian, C. L. Gao, and J.-F. Jia, 
|Phys. Rev. Lett.. 109, 016801 ( 2012). 



[12: 
[13: 

[14 

[is: 

[16 

[17: 
[is: 

[19 

[2o: 

[21 



[22: 

[23: 

[24 



[25: 

[26: 

[27 
[28: 



C.-C. Liu, W. Feng, and Y. Yao, Phys. Rev. LetE] 107, 

076802 (2011). 

B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, 

and A. Kis, Nat. Nanotechnology 6, 147 (2011). 

S. S. Hong, W. Kundhikanjana, J. J. Cha, K. Lai, 

D. Kong, S. Meister, M. A. Kelly Z.-X. Shen, and Y. Cui, 
Nano Lett. 10, 3118 (2010). 

N. Agrait, A. L. Yegati, and J. M. van Ruitenbeek, 

Physics Reports 377, 81 (2003). 

Supplemental Material can be obtained upon request: 

juanjose.palacios@uam.es 

J. M. Krans and J. M. van Ruitenbeek, Physical Review 

B 50, 17659 (1994). 

J. G. Rodrigo, A. Garcia- Martin, J. J. Saenz, and 

S. Vieira, (PhysTRevrLe^ 88, 246801 (2002). 

J. L. Costa-Kramer, N. Garcia, and H. Olin, 

|Phys. Rev. Lett. 78, 4990 (1997). 

R. M. Tahboub, M. E. Guindy, and H. D. Merchant, 

Oxidation of Metals 13, 545 (1979). 

B. J. van Wees, L. P. Kouwenhoven, E. M. M. Willems, 

C. J. P. M. Harmans, J. E. Mooij, H. van Houten, 

C. W. J. Beenakker, J. G. Williamson, and C. T. Foxon, 
Physical Review B 43, 12431 (1991). 

Y. Liu and R. E. Allen, Phys. Rev. B 52, 1566 (1995). 

M. Wada, S. Murakami, F. Freimuth, and G. Bihlmayer, 

Physical Review B 83 121310, (2011). 

M. Mayrhofer-Reinhartshuber, A. Tamtogl, P. Kraus, 

K. H. Rieder, and W. E. Ernst, |J Phys Condens Matter] 

24, 104008 (2012). 

V. D. Das and N. Soundararajan, Physical Review B 35, 

5990 (1987). 

J. J. Palacios, D. Jacob, M. Soriano, A. J. Perez- Jimenez, 

E. S. Fabian, E. Louis, and J. A. Verges, "ALACANT 
quantum transport package", http://alacant.dfa.ua.es. 

D. Jacob and J. J. Palacios, | J. Chem. Phys . 134, 044118 
(2011). 

D. Gosalbez-Martinez, J. J. Palacios, and J. Fernandez- 
Rossier, Physical Review B 83, 115436 (2011). 



